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The motion of metastable helium atoms travelling through a standing light wave is investigated 
with a semi-classical numerical model. The results of a calculation including the velocity dependence 
of the dipole force are compared with those of the commonly used approach, which assumes a 
conservative dipole force. The comparison is made for two atom guiding regimes that can be used 
for the production of nanostructure arrays; a low power regime, where the atoms are focused in a 
standing wave by the dipole force, and a higher power regime, in which the atoms channel along the 
potential minima of the light field. In the low power regime the differences between the two models 
are negligible and both models show that, for lithography purposes, pattern widths of 150 nm can 
be achieved. In the high power channelling regime the conservative force model, predicting 100 nm 
features, is shown to break down. The model that incorporates velocity dependence, resulting in a 
structure size of 40 nm, remains valid, as demonstrated by a comparison with quantum Monte-Carlo 
wavefunction calculations. 
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I. INTRODUCTION 

The dipole force in a standing- wave light field has been 
used to create nanoscale patterns with beams of neu- 
tral atoms for many years. The standing wave can act 
as an array of optical lenses to focus the atoms dur- 
ing deposition onto a substrate, thereby creating a one- 
dimensional structure. The first experiments using this 
technique were performed ten years ago with sodium [J 
and chromium |2| atoms. Since then, atom lithogra- 
phy with aluminium cesium Q| , metastable argon p| , 
and metastable neon p has also been reported. For an 
overview of atom lithography, see Ref . Q . The results of 
a number of these experiments have been compared with 
semi-classical numerical calculations 0, Q . These calcu- 
lations are based on a dipole force that is derived from an 
optical potential to focus the atoms in the standing light 
wave. This force is conservative and does not take into 
account the velocity dependence. Quantum-mechanical 
calculations have been performed as well ^(j- These 
calculations do not suffer from this restriction. How- 
ever, they are computationally intensive and describe the 
atomic motion in one dimension only. 

The purpose of this paper is to investigate in detail 
the 3D atomic motion through a standing-wave light field 
using a semi-classical approach with a non-conservative, 
velocity dependent dipole force. The results of these sim- 
ulations are compared with calculations using the con- 
ventional model with a conservative dipole force. The 
model presented in this paper follows the approach of 
Minogin and Serimaa where the dipole force is rep- 
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resented in the form of a Fourier series, with coefficients 
that are calculated with a continued fraction method. 
Previously, a similar model was used to calculate a one- 
dimensional beam profile of sodium atoms in the far field 
0,0. The model that is presented here describes a full 
3D simulation of the motion of metastable helium atoms 
through a standing light wave for the purpose of atom 
lithography applications. Not only the usual low power 
regime, where the atoms are focused in the standing light 
wave, is investigated, but also a higher power regime. In 
the latter regime the atoms channel through the standing 
wave, undergoing a damped oscillation around the poten- 
tial minima of the light field. This is of special interest 
for lithography, as it allows the production of nanostruc- 
tures with relative insensitivity to the exact alignment of 
atomic beam, substrate, and light field. It is, however, 
also at these higher laser powers that the potential model 
is expected to break down and the inclusion of the ve- 
locity dependence of the force is essential. Experiments 
in this novel regime with metastable helium atoms have 
been performed and, in a future paper, the experimental 
pattern widths will be compared with the calculations 
presented in this paper. With metastable helium atoms, 
nanoscale structures can be created i n a g old film on a sil- 
icon substrate via a two-step process [lj|. First, the high 
internal energy of the focused helium atoms (20 eV) in 
the 2 3 Si metastable state is used to selectively damage 
an organic resist layer through the standing light wave. 
Next, the pattern is transferred to the underlying gold 
film by means of a wet etching process. 

In the simulations, the light field is blue detuned from 
the 2 3 Si -> 2 3 P 2 optical transition (A = 1083 nm) of the 
helium atom. The atoms are therefore attracted to the 
intensity minima of the standing light wave, which min- 
imizes spontaneous emission of photons by the atoms. 
The atomic motion is calculated using only the dipole 
force. The atomic momentum diffusion due to fluctua- 
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tions of the dipole force is neglected. The momentum 
diffusion is caused by the variation in the number of 
absorbed and emitted photons, and by the random di- 
rection of the spontaneously emitted photons. Further- 
more, the atom is assumed to remain in a steady state, 
i.e., transient effects of the dipole force are neglected. 
In order to investigate the effects of momentum diffu- 
sion and transient effects, some results from Monte Carlo 
Wave Function (MCWF) calculations similar to Lee [ljj 
are presented as well. These calculations include the 
velocity dependence of the force, momentum diffusion, 
atomic diffraction as well as transient effects. However, 
they are one-dimensional and require significant compu- 
tational resources. 

In the next section, a general expression for the dipole 
force is derived, and both semi-classical models are out- 
lined. Sec. IIIII describes the method and the parameters 
of the numerical simulation. The results of the semi- 
classical simulations are presented in Sec. II VI and they 
are compared with the MCWF simulations. Finally, con- 
cluding remarks are given in Sec. Ivl 



II. THEORETICAL MODELS 
A. Dipole force of a standing wave 

The light force experienced by a two-level atom due to 
the presence of a classically described light field can be 
determined by calculating the change in momentum p of 
the atom under influence of the Hamiltonian: 



H = H A -d- E, 



(1) 



where Ha contains the internal and kinetic energy of the 
atom and d ■ E is the dipole interaction operator that 
describes the coupling between the atomic dipole moment 
d and the electric field component E of the light field. 
According to the Ehrenfest theorem, the quantum-me- 
chanical analogue of the radiation force is given by the 
Heisenberg equation of motion 



dip) 
dt 



V(d-E) 



VE. 



(2) 

In the last step of Eq. the expectation value of the 
electric field operator E is replaced by the value at the 
atomic centre of mass. This is legitimate in the electric 
dipole approximation, where the wavelength of the light 
field A is large compared to the de Broglie wavelength 
Ads = fo/\p\ of the atom, and spatial variations of the 
electric field on the scale of the atomic wave-packet can 
be neglected. 

The expectation value of the electric dipole operator 
d can be written in terms of the atomic density matrix 
p, which describes the quantum-mechanical state of the 



two-level atom, as 

(d) = Tr(pd) = d ge (p ge + p eg ) 

— 2d ge (u(t) cos tut — v(t) sin cut), 



(3) 



where the atomic density matrix elements p ge and p eg — 
p* e are the electronic coherences between the ground 
state and the excited state of the atom, oj is the fre- 
quency of the radiation field, and u(t) and v(t) are two 
components of the Bloch vector. In the rotating-wave ap- 
proximation, where non-resonant terms of the atom-light 
interactions are neglected, the components of the Bloch 
vector can be written as 



u{t) 



«(*) = iliP: 



+ Pege ), 



w(t) 



2 (Pee — Pgg)i 



Pegt 



(4) 



where the atomic density matrix elements p gg and p ee 
are the populations of the ground state and the excited 
state of the atom, normalized to p gg + p ee = 1. For a 
travelling wave, the electric field component of the light 
field is given by 



E(r, t) = e(r) E (r) cos {tot + $(f)), 



(5) 



where e*(r), Eo(r), and (= —k ■ r) are the polariza- 
tion, amplitude and phase of the light wave respectively 
at the atomic centre-of-mass position r = (x,y,z). The 
Rabi frequency f2(r) of the light field is then defined as 



d ge ■ e(r) E (f) 



(6) 



The general expression for the light force (in the electric 
dipole and rotating-wave approximations) can be written 
as [lU 

?) = F{r) = -hu st VSl(f) - hn(r)v st V<P(r). (7) 

The two parts on the right-hand side of Eq. (jJJ) are the 
dipole force, proportional to the gradient of the Rabi fre- 
quency Cl(r), and the scattering force, proportional to 
the gradient of the phase <&(r) of the light field. The 
Bloch vector components u(t) and v(t) are replaced by 
their time-independent steady-state values u s t and v s t re- 
spectively. This is valid in the adiabatic approximation, 
where the atom moves slowly enough in the light field to 
maintain an equilibrium between its internal state and 
the radiation field. 

The steady-state values can be found by solving the 
equations of motion of the optical Bloch vector that de- 
scribe the time evolution of a two-level atom in a light 
field 
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(8) 
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where V /2ir is the natural linewidth of the excited state 
and A/27T is the detuning between the light field fre- 
quency lo and the frequency of the atomic transition. 

The electric field component of a standing-wave light 
field, composed of an incident and back-reflected Gaus- 
sian wave travelling in the ^-direction, can be written as 



Eq. H14fl is a well-known expression for the dipole force 
and it is commonly used for semi-classical calculations 
of atomic motion in a standing- wave light field 0, 0, E| . 
Since this force is conservative, the kinetic energy of the 
atom at any moment is determined by the local potential 
of the light field. 



E(r, t) = eEo cos (wt) sin (kx) exp 



(9) 



where k is the wave number of the light and w y and w z 
are the waists of the Gaussian beam profile in the y and 
z direction respectively. Since this electric field has no 
phase dependence, the phase gradient term in Eq. (J7J) 
vanishes. The general solution of the light force in a 
standing-wave light field then becomes 



F sw (r) = -hu st VQ(r). 



(10) 



The steady-state value of the optical Bloch vector com- 
ponent u s t can be derived using two different approaches; 
one that neglects the atomic velocity f Sec. Ill B"|) . and one 
that is valid for arbitrary atomic velocities (Sec. HI Cjl . 



B. Atom at rest — potential model 

When the transverse velocity of the atom in a standing- 
wave light field is negligible (k ■ v -C T), the atom travels 
over a very small distance compared to the optical wave- 
length A during the relaxation time T _1 of the atom. The 
optical Bloch equations given in Eq. (JSJ can then be con- 
sidered as a set of coupled linear differential equations 
with time-independent coefficients. The steady-state so- 
lutions are found analytically by setting u = it = w = 
and they are given by 



where 
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is the saturation parameter. The final expression for the 
dipole force acting on an atom at rest in a standing-wave 
light field, can now be found by combining Eq. 110(1 and 
Eq. (TU 



-hA 



20 2 (r)+r 2 + 4A 2 ' 



(13) 



This force is conservative and it can be written as the 
gradient of a potential (IH H(| 



F pot (f) = -VU(r) 



-V^ln[l + ,s(f)]. (14) 



C. Moving atom — Minogin model 

When an atom moves with velocity v in the 
standing-wave light field, the position of the atom f be- 
comes explicitly time dependent as r(t) = vt. However, 
the time dependence of the transverse coordinates y and 
z can be neglected, since the wavelength of the light A 
is much smaller than the waist of the Gaussian beam 
profile. This means that the optical Bloch vector com- 
ponents u, v, and w change more rapidly along the axis 
of the light field (^-direction) than in the transverse di- 
rections y and z. Therefore, for finding the steady-state 
solution of the optical Bloch vector components, only the 
time dependence of the ^-coordinate has to be taken into 
account. The Rabi frequency 0(f) can then be written 
as a periodic function of time 



n(f,t) = Clo(y, z) sin(fcx) 

= n (y,z)sin(kv w t), 



(15) 



where Qo(y,z) is the peak Rabi frequency at the anti- 
nodes of the standing wave, and v x is the velocity of the 
atom parallel to the axis of the light field. With the Rabi 
frequency given by Eq. I(15J1 . the coefficients of the cou- 
pled linear differential equations in Eq. (JSJ become time 
dependent. Consequently, the optical Bloch equations 
can no longer be solved analytically. However, a steady- 
state solution of u, v and w can be found by expanding 
each of them in a Fourier series 



(16) 



where the common notation h is used for u, v, and w, and 
the quantities h n satisfy the reality condition /i_„ = h* n . 
By substituting the Fourier expansions Eq. 1(16(1 and the 
expression for the Rabi frequency given by Eq. (II 51) in 
the optical Bloch equations Eq. (JSJ , a set of recursive al- 
gebraic equations is obtained for the optical Bloch vector 
components ^l| 



(r/2 + inkv x )u r , 
(r/2 + inkv x )v„ 



-fi (W n -i + W n+ i) - Au n , 



(T + inkv x )w n = Q (u„_i + v n+ i) - — 5 n0 . (17) 

The dipole force is now also represented in the form of a 
Fourier series, where the Fourier coefficients can be cal- 
culated from the optical Bloch vector components. The 
expression for the dipole force can be separated into com- 
ponents parallel and perpendicular to the axis of the 
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standing-wave light beam (x-direction) as 

F\i(r) = F°(r) (18) 

OO 

+ X^ 1 )" (^u(r) cos2nkx + F^ l (f) sin2nfcx) , 

n=l 

and 

F ± (r) = F°(r) (19) 

OO 

+ £(-l) n (jxS(» s ) cos2nkx + Fl n v {r) sin2nfca;), 
n=l 

where the coefficients of the Fourier series are given by 
F^(f) = -2hkfl lmui, 

^|f™(0 = -2hkn Q lm(u 2 n+l ~ «2n-l)i 

fjf o n (»0 = -2ftfcO Re(u2„+i -u 2n -i), (20) 

and 

^(r) = 4h(^r + A; J HoReui, 
Fft (f) = 4fi f + 4^ Re (u 2n+1 + u 2n ^), 

The coefficients u n are obtained from the solution of the 
recursion relations (Eq. I|17[ll in the form of convergent 
continued fractions. The dipole force given by Eq. i|18|) 
and Eq. (|TT)|l is dissipative, and the atomic energy is 
therefore not conserved. The momentum change of the 
atom can be attributed to an additional damping or heat- 
ing force, dependent on the intensity of the light field and 
on the sign of its detuning from atomic resonance. 

III. NUMERICAL SIMULATIONS 
A. Calculation method 



described in detail by Minogin and Letokhov [17| and it 
is summarized for completeness in the appendix. The so- 
lutions for the coefficients are found in the form of con- 
vergent continued fractions. For the simulations, these 
continued fractions are calculated with an accuracy bet- 
ter than 10 -16 , which requires a maximum of 5000 terms. 
This precision is necessary to calculate the coefficients u n 
with an accuracy better than 1 ppm. The Fourier series 
that represent the expression for the dipole force contain 
at most 2500 Fourier coefficients. A convergence analy- 
sis has shown that the dipole force change is negligible 
at higher accuracies of the coefficients u n and when more 
terms are added to the Fourier series. Therefore the level 
of convergence is assumed to be sufficient. 



B. Parameters 

The calculations are performed on a beam of helium 
atoms in the 2 3 Si metastable state, which has a lifetime 
of about 8000 s. The atomic beam has a mean longi- 
tudinal velocity of 2000 m/s and a longitudinal velocity 
spread (full width at 1/e 2 height) of 650 m/s. After col- 
limation of the beam, the transverse velocity spread of 
the atoms is reduced to about 3 m/s. For the calcula- 
tions, each atom is assigned a longitudinal and transverse 
initial velocity that is randomly picked from Gaussian 
velocity distributions with the above described averages 
and spreads. 

The light of the standing wave has a wavelength of 
1083 nm, driving the 2 3 Si — > 2 3 P2 optical transition of 
the helium atom, which has a natural linewidth r/27r = 
1.6 MHz. By detuning the light field relatively far to the 
blue side of the atomic resonance (A/2ir — 375 MHz), 
the atoms are attracted to the nodes of the standing- 
wave light field, which reduces spontaneous emissions. 
An upper limit for the detuning forms the 2 3 Pi state, 
which energy level lies 2.3 GHz above the 2 3 P2 state. 
For very large blue detunings of the light field from the 
2 3 P2 state, the atom can thus interact with the light 
field via the 2 3 Si -> 2 3 Pi transition. 

The Rabi frequency can be calculated from the inten- 
sity of the light field as 



The dipole force derived from the two models de- 
scribed in Sec. Ill Bl and Sec. Ill Cl is used to calculate the 
atomic movement through a standing- wave light field. By 
straightforward numerical integration of the Newtonian 
equations of motion, the change of the atomic velocity 
and position under influence of the dipole force at the 
current position of the atom in the light field is calcu- 
lated. In this way a full 3D simulation of the atom tra- 
jectories and velocities in the standing-wave light field 
is performed. The atomic pattern formation is mapped 
with 2D atomic distribution plots in the xy-plane and in 
histograms. 

The calculation of the coefficients u n , required for the 
calculation of the dipole force in the Minogin model, is 



Q(?\-r* (22) 

where I sat = 0.17 mW/cm 2 is the saturation intensity of 
the optical transition and J(r) is the intensity profile of 
the standing- wave light field, given by 

I{r) = I sin 2 (kx) exp (~% - (23) 

The Gaussian light beam has a circular beam profile with 
a waist (1/e 2 radius) w y — w z — 331 /im. The quantity 
Iq is the intensity of the light field at the anti-nodes of 
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Figure 1: Velocity dependence of the re-component of the dipole force of a standing- wave light field in the focusing regime (left) 
and the channelling regime (right) at position (a;, y, z) = (3A/8, 0, 0) in the standing wave and at a detuning of A/2n = 375 MHz. 
The dipole force according the Minogin model (solid lines) is normalized to the dipole force from the potential model (dashed 
lines) . 



the standing wave and it is given by 



h = 



SPp 

■KWyW z ' 



(24) 



where Pa is the power of the incident light beam. De- 
pending on this power, two different regimes can be dis- 
tinguished for guiding the atoms through the standing- 
wave light field. At low power, the atoms can be focused 
at the centre of the Gaussian light beam. For large de- 
tunings (A 3> r), the power required for this focusing 
can be calculated from 



Pn = 5.37 



2hT 2 k 2 ■ 



(25) 



This focusing power is independent of the waist of the 
light beam. For the conditions mentioned above, Pq = 
2.4 mW. For high-power light fields, the atoms oscillate 
through a potential minimum of the standing wave and 
the sign of the transverse velocity of the atoms changes 
many times. This is called the channelling regime. For 
calculations in this regime a power Pq = 800 mW is 
used. The dipole force in the channelling regime is then 
about one order of magnitude larger than in the focusing 
regime. 



C. The dipole force 

The velocity dependence of the ir-component of the 
dipole force in the focusing and channelling regime for 
both models is depicted in Fig.^ The graphs show that, 
in the focusing regime (left graph), the difference be- 
tween the dipole force according to the potential model 
and the Minogin model is very small (at most 0.3%). 
In the channelling regime (right graph), however, the 



dipole force of the two models differs significantly. The 
force shows a negative slope for small transverse velocities 
(|fcu x /r| < 0.2), which means that the force is a damping 
force. This cooling effect is essentially the "blue-detuned 
Sisyphus cooling" introduced by Dalibard and Cohen- 
Tannoudji 18] . At larger velocities, the dipole force in 
the Minogin model is significantly larger than in the po- 
tential model. It is therefore expected that the distinc- 
tion between the two models will be most pronounced in 
the channelling regime. Furthermore, in the channelling 
regime, the Minogin model shows some resonance peaks 
at high atom velocities. These are called Doppleron res- 
onances |19| . and occur when the atom is excited to the 
2 3 P2 state by multiple photon absorptions and emis- 
sions. When an atom absorbs (n + 1) photons from one 
wave of the standing-wave light field, and emits n pho- 
tons into the other, these resonances appear at velocities 



kv x /T = ± 



A 



(2n+ l)r' 



(26) 



IV. RESULTS 

Fig.[2]shows a projection on the xz-plane of the atomic 
motion through the standing-wave light field in the focus- 
ing regime (upper graphs) and channelling regime (lower 
graphs). The centre of the Gaussian light beam, which 
propagates in the rr-direction, is located at z = 0. The 
graphs cover half a wavelength (A/2 = 542 nm) of the 
standing-wave light field, showing exactly one potential 
minimum through which the atoms travel in the positive 
^-direction (from left to right). To allow the atoms to ex- 
perience a maximum dipole force, the starting y-position 
is y — for all atoms. However, since the simulations 
are 3D and the atoms have an initial velocity in the y- 
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Figure 2: Atom trajectories through the standing-wave light field in the focusing regime (upper graphs) and the channelling 
regime (lower graphs). Semi-classical calculations are performed with the potential model (left) and the Minogin model (centre). 
The right graphs show atomic density distributions of Monte Carlo Wave Function (MCWF) simulations. 



direction, this position is not maintained while travelling 
through the light field. 

Since the velocity dependence of the dipole force in 
the focusing regime is almost negligible (see Fig.^), the 
atomic motion in this regime according to the potential 
model (upper left graph) is very similar to the calcula- 
tion with the Minogin model (upper centre graph). In 
the channelling regime, the calculated trajectories differ 
considerably for the potential model (lower left graph) 
and the Minogin model (lower centre graph), due to the 
higher average value and the strong velocity dependence 
of the dipole force in the latter model (see Fig. [TJ . The 
potential model is invalid in this regime, since the atoms 
have a transverse velocity spread of 3 m/s (kv x /T = 1.7), 
and the transverse velocity of the atoms can even be- 
come on the order of 10 m/s in the standing wave. In 
the channelling graph of Fig. |3 the additional damping 
force results in an increased confinement of the atoms in 
the potential well of the standing wave. The upper right 
and lower right graphs of Fig. [21 show one-dimensional 
MCWF calculations for the focusing regime and chan- 
nelling regime respectively. In these calculations atomic 
diffraction, velocity dependence and momentum diffusion 
are fully incorporated. The plots show the atomic den- 
sity distribution on a grey-scale for 225 quantum tra- 
jectories, corresponding to the same initial distribution 
as used for the semi-classical simulations. The charac- 
teristic feather-like structures in the channelling regime 
(lower right graph) do not constitute a specific quantum 
feature. They correspond to (approximate) caustics in 
the semi-classical trajectories and they become clearly 
visible if a larger number of trajectories is plotted in the 
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Figure 3: Distribution plots of the positions where the atoms 
hit the sample after travelling through the standing-wave light 
field at sample position z — w z /<i = 83 fim. From left to right, 
the first to graphs show plots of the focusing regime using the 
potential model (first graph) and the Minogin model (second 
graph). The last two graphs show the channelling regime us- 
ing the potential model (third graph) and the Minogin model 
(fourth graph). 



same way as in the lower centre graph. 

Distribution plots of the atoms are shown in Fig. [3] 
Every dot in these graphs represents a position where 
the atom has hit the sample after travelling through the 
standing-wave light field. The calculations are performed 
on atoms in a lattice of 101 x 101 atoms. The starting 
positions of the atoms range from (x,y) = (—A/2, —2w y ) 
to (x,y) — (A/2, 2w y ). The plots only show half of this 
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Figure 4: Histograms of the atom distributions taken at z = w z /4 = 83 /im. The upper graphs show the distributions for 
the focusing regime in the potential model (left), Minogin model (centre), and the MCWF simulations (right). Similarly, the 
channelling regime is represented by the lower graphs. 



window in the ^-direction, from x = —A/4 to x = A/4. 
In the focusing regime, the light-field does not act as a 
perfect lens, but shows some abberations. This imper- 
fection is enhanced by the longitudinal and transverse 
velocity spread of the atoms. As a result, the best posi- 
tion to place the sample is not at the centre of the light 
beam (z = 0), but slightly behind it. For the channelling 
regime, the sample position is far less critical, but the 
best results are at similar positions. Therefore, the com- 
parison of the two regimes for the different models in 
Fig- El is performed at z = w z /A = 83 /im, where optimal 
results are expected. 

The first two plots show the atom distributions in the 
focusing regime for the potential model (first plot) and 
the Minogin model (second plot). As expected from the 
comparison of the trajectories in this regime, shown in 
the upper graphs of Fig. |3 the differences between the 
two models are negligible. The distributions in the chan- 
nelling regime with the potential model (third plot) and 
the Minogin model (fourth plot) show that the atoms are 
confined for a wider range along the y-axis as compared 
to the focusing regime. Furthermore, for the Minogin 
model, the atoms are more localized to the nodes of the 
standing wave. 

To make a more quantitative comparison of the four 
plots of Fig. histograms of the atom distributions are 
presented in Fig. Only the atoms located between 
y = —100 /im and y — 100 /im are taken into account 
for the histograms, because atoms outside this region 
do not contribute to the desired pattern (in the focus- 
ing regime). The black curves through the histograms 
are Lorentzian fits from which the Full Width at Half 



Maximum (FWHM) of the distribution can be deduced. 
The upper graphs, that represent the calculations in the 
focusing regime, show that the FWHM of the distribu- 
tion is 150 nm for both semi-classical models as well as 
for the quantum-mechanical model. For the calculations 
in the channelling regime (lower graphs), the distribu- 
tion of the potential model (left graph) has a FWHM of 
98 nm. This width deviates clearly from the distribu- 
tion of the Minogin model (centre graph), which has a 
FWHM of 40 nm. The smaller width of this distribu- 
tion can be explained from the additional cooling force 
that is included in the Minogin model. The distributions 
for the full quantum-mechanical calculations are shown 
in the right graphs. The FWHM of the central peak is 
45 nm, somewhat larger than the results from the con- 
tinued fraction semi-classical results. This larger width 
is mostly due to the fact that the steady-state situation 
for the atoms is not completely realized during the inter- 
action time. Achieving a steady state requires a number 
of spontaneous emissions per atom. Due to the large de- 
tuning of the light field from atomic resonance, and due 
to the fact that the atoms are channelled in the region 
of low light intensity, almost one half of the atoms never 
undergoes a spontaneous emission. The influence of dif- 
fusion is seen in the wings of the MCWF distributions. 
However, overall the MCWF results agree very well with 
the semi-classical calculations using the Minogin model, 
confirming the large reduction in the FWHM of the dis- 
tribution due to the cooling force. 
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V. CONCLUSIONS 

The simulations have shown that nanoscale patterns 
can be created utilizing the optical dipole force to guide 
atoms through a standing-wave light field. For the cal- 
culation of the atomic motion, two models have been 
applied: a simple one that uses the conventional dipole 
force, derived from a potential, and another one that in- 
cludes the velocity dependence of this dipole force. In the 
'conventional' focusing regime, where a low-power light 
field is used, the differences between the potential model 
and the Minogin model arc negligible. Both models show 
that in principle nanoscale pattern sizes of 150 nm can be 
achieved with the specified experimental parameters. It 
should be noted, that these parameters are not optimized 
for focusing: using a smaller laser focus and more laser 
power, much tighter focusing can be achieved without en- 
tering the channelling regime. However, the sample po- 
sition is very critical, and good focusing is only achieved 
for atoms that pass the light field close to its centre. The 
high-power channelling regime, that we are presently us- 
ing for nanostructure production with metastable helium 
atoms, is more robust. It is therefore better suitable for 
creating narrow patterns. In this regime, the Minogin 



model shows strong deviations from the potential model, 
due to the contribution of velocity dependent terms to 
the dipole force. Because of the large transverse veloc- 
ity (kv x > r) of the atoms in the standing-wave light 
field, the potential model breaks down and should not 
be applied for calculations in this regime. The results 
of the calculations with the Minogin model show struc- 
tures with a FWHM of 40 nm. This model does not 
take into account the momentum diffusion and transient 
effects, which leads to broadening of the structure size, 
as is shown by the MCWF simulations. However, the 
MCWF simulations support the results of the Minogin 
model by showing a better confinement of the atoms in 
the standing-wave light field than expected from calcu- 
lations with the conventional model. This is due to an 
additional cooling force that is not incorporated in the 
conventional model. 
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Appendix A: CONTINUED FRACTION METHOD 

The determination of the dipole force in the Minogin 
model, given by Eqs. (|18l) and l|19l) . requires the opti- 
cal Bloch vector components u n , that are embedded in 
the expressions of the Fourier coefficients in Eqs. I|20|) and 
(|21|l . The coefficients u n can be calculated with a contin- 
ued fraction method according to Minogin and Letokhov 
[l7| . A summary of this method is described below. 

The upper expression of Eq. I|17|) can be rewritten in 
the form 



r/2 + inkv x n 



(Al) 



Since v n is non-zero for odd n, and w n is non-zero for 
even n, the middle and lower expressions of Eq. (|17|) can 



be combined to a single expression as 



b n - D n (b n -i + b n+ i) 



where 



K = 



w r , 



for odd n, 
for even n, 



(A2) 



(A3) 
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and the coefficients D n are given by 

!fl (r/2 + inkv r ) r , , 
' ■> ; „ for odd n, 
{T/2 + inkv x ) 2 + A 2 
-n i ^ i — for even n. 
1 + inkv x 

The coefficients 6„ fulfil the reality condition 



b -n = b n , 



(A4) 



(A5) 



and therefore only non-negative integers (n > 0) have 
to be considered. The relation between two successive 
quantities b n and 6„+i can then be found by substituting 



bn+l = qnbn 



(A6) 



into Eq. (| A2|) . This leads to a recursion relation for q n 
given by 



1 



1 



In 



(A7) 



The quantity qo can now be expressed as a convergent 



continued fraction 



qo = 



Pi 



P2 



P3 



where 



The quantity bo 
n = and gives 



W 



bo 



i 



4ft/r Re (g ) - 2' 



(A8) 



(A9) 



is found by solving Eq. (| A2|) for 



(AlO) 



From Eq. I|A6|) and this initial coefficient bo, all succes- 
sive coefficients b n can be calculated. By substituting 
the coefficients b n (= v n for odd n) into Eq. (| Al|) . the 
coefficients it n can be found straightforwardly. 



